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Within the relativistic mean-field approach, we extend the Miyazaki model, where the 
C^l . NNcr and NNoj interactions are modified to suppress the couplings between positive- and 

» I ■ negative-energy states of a nucleon in matter. Assuming appropriate density-dependence 

of the meson-nucleon couplings, we study nuclear matter and finite nuclei. The model can 
reproduce the observed properties of 16 O and 40 Ca well. We also examine if the model is 
natural. 

Recently, the relativistic mean-field approach with density-dependent meson- 
nucleon couplings draws much attention.^ It is an effective model for the Dirac- 
Brueckner-Hartree-Fock (DBHF) theory,^ which can reproduce the saturation prop- 
erty of nuclear matter using the one-boson exchange potentials extracted from the 
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nucleon-nucleon scattering data. In the DBHF calculation, the relativistic effect 
provides a strong density-dependent repulsion, which is originated from the nucleon- 
antinucleon pair term (Z graph), and it is vital to obtain the nuclear saturation 
property. It should be noticed that a nuclear model based on the quark substructure 
of a nucleon, for example, the quark-meson coupling (QMC) model,® the quark-mean 
field (QMF) model, 4 ' also gives density-dependent meson-nucleon couplings through 
the scalar field in a nuclear medium, namely the scalar polarizability. 3 ' Thus, it seems 
quite natural that the meson-nucleon couplings depend on the nuclear environment. 

About a decade ago, Miyazaki has proposed an interesting, relativistic mean- 
field model for nuclear matter, in which the NNcr and NNw vertices are modified to 
reduce the couplings between positive- and negative-energy states of the in-medium 
nucleon (the H — couplings). Although the H — couplings play an important role in 
the relativistic nuclear models including nucleon-nucleus (NA) scattering (with the 
relativistic impulse approximation (RIA)) at intermediate energies, it is known that 
the effect of the coupling to negative states is too strong to produce the NA scatter- 
ing observables at low energies P 1 Tjon and Wallace have remedied this problem by 
developing a generalized RIA, in which the different +— couplings from the usual 
RIA are introduced.® The vertex modification studied by Miyazaki^ may enable us 
to include such variation of the H — couplings at the relativistic mean-field level. The 
modified vertices finally result in the density-dependent NNcr and NNu; couplings, 
which can simultaneously reproduce the nuclear matter properties and the Dirac 
scalar and vector optical potentials given by the DBHF calculation. 

In this Letter, we generalize the Miyazaki model, and study not only the nuclear 
matter properties but also single-particle energies of finite nuclei. Lastly, we discuss 
naturalness of the modelP 



typeset using pTpT^.ds (Ver.0.9) 
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We now modify the vertices of NNu and NNcj couplings using the energy projec- 
tion operators, A^(p) = + M)/2M , where p is the four-momentum of a nucleon 
and M is the mass. Since the vertex, r(= I or 7^), is expressed by 

r = A + {p')rA + ( P ) + A-{p')rA~( P ) + A + {p')rA~( P ) + A~( P ')rA + ( P ), (i) 

it may be possible to vary the strength of the +— couplings, introducing two param- 
eters, < Ai,A 2 < 1, asP 

r -» \ 1 [A + { P ')rA+( P ) + A-{p')rA-{ P )} + \ 2 [A + {p')rA-(p) + A-( P ')rA+(j>)},{2) 
+ (\ 1 + x 2 )rM 2 



2M 2 



(3) 



In the original Miyazaki modelP Ai is chosen to be unity for the scalar (I) 
vertex, while A2 is unity for the vector (7^) vertex, because the parameters are 
supposed to be constants. However, in general, the strength of the +— couplings 
may depend on the nuclear environment through the Pauli blocking, Z graphs etc. 8 ' 
To take account of those effects in the model, we here suppose that A simply depends 
on the nuclear density, p v : 

* = 1 -(£)'' (4> 

where po is the saturation density, and each A has two parameters, a and b. Note 
that, in the limit p v — > 0, r is identical to the original form Eq.([T]). 

Using the vertex (J3j) and the mean-field approximation for the meson fields, the 
Lagrangian density is given byS 

£ = M)ip - \™lo 2 + \mloj 2 

+ [(A! - Ag)$? W) + (Af + \ s 2 )M 2 ^]a 

■ Lj w-a^?) 7 ^) + (a? + a^mW#^ (5) 
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where a and u are respectively the mean-field values of the a and oj mesons, and 
A| (i = 1,2) is the parameter for the scalar (vector) vertex. The meson mass 
and the NN<j(o;) coupling constant in vacuum are respectively denoted by m ff y and 
9<t(u>) ■ 

Following the prescription explained in RefP , we can construct an effective 
Lagrangian density, in which the effect of variation of the +— couplings in matter is 
included, 

C e & = ~ 1°U V - M*)^ - l -m 2 y + \mluj 2 , (6) 

where the effective nucleon mass, M*, and the Dirac scalar, U s , and vector, U v , 
potentials in matter are defined as 

M-M*= g* a a = -U s , U v = g*u, (7) 
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with the effective coupling constants 
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[(Af + AS) + (Af 

(A? 



2 

~[(Af + A?0 



A|)(m* 2 -?; 2 )] 5(7 , (8) 

A^m* 2 -*, 2 )]^, (9) 
1, the effective coupling 
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m* = M*/M and v = U v /M. Note that, when Af 
constants coincide with those in the Miyazaki model. 

The energy per nucleon, W, for symmetric nuclear matter is then written by 

1 2 



W = ^E* F + \m*^ + U v 
4 4 p v 



2M 



M + 



2M 

c7p 
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A? + A| + (Af 

2 



AI)(m 



* 2 



Af + A^ 



(X\ - X v 2 ){m* 2 - v 2 ) 



(10) 



where E* F = (k F + M* 2 ) 1 / 2 (fcp the Fermi momentum) 

Q . . — „2 „_ /^,2 A/f n„rl ^ — f /\/T* /m-2"\ fl._ I?* 7\/f*- 



M«) = ^(^Po/m^^M, and /? s = (AfVvr 2 )[A: F ^ 
scalar density). 

From the self-consistency conditions, (3VF/9m* 
the meson fields should satisfy, one finds 

4c 
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p v = 2k F /37r 2 , p = p v /po, 
M* 2 \n{(k F + E* F )/M*)\ (the 



and (dW/dv) = 0, which 
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where 

a s = Af + A| + (Af 
Qo> = Af + A2 - (Af 



*2 



Af + A^ 



A5)(m 
\l)(m* 2 -v 2 ), 

-X v 2 )(m* 2 + v 2 ) 



AI 



-(Af 
(Af-A|)(m* 2 



2(Af - \ v 2 )m*v, 



2m* + v 2 ) + 2(Af 



c=[Af + A| - (Af 
+ 4(Af-A|)(Af - 



-A|)(m* 2 -2m* + t> 2 )][Af + A'^ 
X v 2 )(l-m*)m*v 2 . 



AI)(1 
- (Af 



m }v- 



(11) 

(12) 
(13) 
(14) 

(15) 

I 

(16) 



Giving the A parameters, the coupling constants in vacuum, g a and g^, are 
determined so as to fulfill the saturation condition, dW / dp\p=\ = and W(p = 1) = 
— 15.75MeV (po = 0.15 fm~ 3 ). Using those coupling constants, one can calculate the 
effective nucleon mass (m*) and the vector potential (v) at any density. 



In the present calculation, we, however, set that Af 



Ao 



1 



aA(pv/po) bA 



and A| = Af 



1 — ca(Pv/ Po) dA to reduce the total number of parameters. (We 
call this "type A".) Thus, for type A the density-dependence of the NNu coupling 
is identical to that of the NNu; coupling, i.e., g* a jg a = gZ/dui (see Eqs.dH) and ([9])). 
This choice may be justified, because it has already been found in Ref.^ that the 
case of Af = X 2 and X 2 = Af gives the best result for the nuclear matter properties. 
In contrast, we shall also study an alternative: Af = Af = 1 — cib(Pv/ Po) bB and 
X 2 = X 2 = 1 — cb(Pv/ Po) ds ■ (We call this "type B".) Thus, each type eventually has 
four parameters (a ~ d) to fit the observed data. 
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Table I. Parameter sets for type A and B. The 
nuclear incompressibility K (in MeV) is 
also shown. 



Now we are in a position to show 
our numerical results. We determine the 
coupling constants, g a and g^, so as to 
exactly reproduce the saturation condi- 
tion of nuclear matter. In addition, the 
parameters, a ~ c£, are tuned so as to 
produce the scalar U s and vector U v po- 
tentials of the DBHF calculation and the observed incompressibility (K = 210 ± 30 
MeV) as precisely as possible. 

We find that g a = 12.24(16.97) and g w = 14.97(20.69) for type A (B). In Table[H 
we list the parameter sets for type A and B. In Figs. 1 and 2, the scalar and vector 
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d 
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0.1 
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0.40 


0.3 


0.37 


0.3 


202.9 
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Fig. 1. The scalar potential. The dashed, solid 
and dotted curves are, respectively, for type 
A, B and QHD. The DBHF result is shown 
by solid squares. 




P, 7 Po 



Fig. 2. The vector potential. The curves are 
labeled as in Fig. 1. 



potentials calculated with the present parameter sets are shown, together with the 
results of DBHF calculation 2 ' and Quantum Hadrodynamics (QHD)P The DBHF 
result is well reproduced by the present model up to p v /po = 2.0 ~ 2.5. We again see 
from TableUthat g^/da — dZ/d^ even in type B because bs = ds and ae ~ eg. Thus, 
for fitting the scalar and vector potentials of the DBHF calculation and the observed 
incompressibility simultaneously, it may be favorable that the density-dependence of 
the NNct interaction is very close to that of the NNu interaction. 



Table II. Binding energy per nucleon W (in MeV), rms charge radius r c h (in fm) and difference 
between nuclear radii for neutrons and protons r n — r p (in fm). The QHD result is also included. 
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463.4 


3.482 


5.78 


-0.080 


2.83 


4.22 


-0.047 
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3.482 


7.33 


-0.074 


2.77 


6.06 


-0.041 


QHD 


523.8 


3.482 


6.24 


-0.055 


2.75 


4.85 


-0.033 


Exp. 




3.482 


8.45 


0.05±0.05 


2.73 


7.98 
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Table III. Model predictions for the energy spectrum of Ca. 
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Fig. 3. Charge density distribution for 40 Ca compared 
with that of QHD. The curves are labeled as in Fig. 1. 



For finite nuclei, Eq.© 
gives a set of coupled non-linear 
differential equations, which 
may be solved by a standard 
iteration procedure.^ For ex- 
ample, we have calculated the 
properties of 16 O and 40 Ca, and 
the result is presented in Ta- 
ble [H] In the calculation, we 
have adjusted the a mass (m CT ) 
so as to yield the observed root- 
mean-square (rms) charge ra- 
dius of 40 Ca: r ch ( 40 Ca) = 3.48 
fm. In Table IIIIl we give 
the single-particle energies for 
40 Ca. We can see from the table 
that the model, especially type 
B, produces the good result. The spin-orbit force in the present model is thus suffi- 
cient to reproduce the observed energy levels and it is comparable to that of QHD. 
The charge density distribution for 40 Ca is also illustrated in Fig. El together with 
the QHD result. Note that the observed distribution, which is not shown here, is 
very close to that of type B. 

Finally, we shall examine the present model using Georgi's "naive dimensional 
analysis" (NDA).^'^ i n general, an effective field theory at low energy will con- 
tain an infinite number of interaction terms, which incorporate the compositeness of 
the low-energy degrees of freedom, i.e., hadrons, and it is then expected to involve 
numerous couplings which may be nonrenormalizable. The NDA gives a systematic 
way to manage such complicated, effective field theories. After extracting the di- 
mensional factors and some appropriate counting factors using NDA, the remaining 
dimensionless coefficients are all assumed to be of order unity. This is the so-called 
naturalness assumption. If theory is natural, one can then control the effective La- 
grangian, at least at the tree level. In the present case, the model involves the NNcr 
and NNw interactions. Using NDA, we then find that the dimensionless coefficients 
corresponding to those couplings are all smaller than 2.0. Thus, the model is natural. 
In summary, we have extended the Miyazaki model, - where the NNcr and NNw 
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couplings are modified to suppress the H — couplings, and studied two (A and B) 
types of density-dependent meson- nucleon vertices. Assuming an appropriate density 
form at the vertex, the parameters are adjusted so as to produce the scalar and 
vector potentials of the DBHF calculation and the observed incompressibility K 
as precisely as possible. The density-dependence such as g%/g<j ~ gZ/du ma Y be 
eventually favorable for fitting the DBHF result and the observed nuclear data. 
Using such coupling constants, we have studied the properties of nuclear matter and 
finite nuclei ( 16 and 40 Ca). The model can reproduce the experimental data well. 
Furthermore, NDA tells us that the model is natural. It is thus vital to include 
appropriate density-dependence of the meson-nucleon interactions in the relativistic 
mean-field approach, which may be attributed to many-body effects (like the Pauli 
exclusion) and the quark substructure of an in-medium nucleonP 
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